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Abstract 

T— I , 

f*"*) . We show that the coined quantum walk on a line can be understood 

£N| ' as an interference phenomenon, can be classically implemented, and in- 

deed already has been. The walk is essentially two independent walks 

■ associated with the different coin sides, coupled only at initiation. There 

is a simple analogy between the evolution of walker positions and the 
propagation of light in a dispersive optical fiber. 
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The quantum random walk (QW) was first proposed ten years ago by Aharonov 
et al. pQ as the quantum analog of the classical random walk (RW) . QWs are 
receiving much attention as some problems are best solved in classical 

computation with algorithms based on RWs, it is expected that these type of 
problems could be solved even faster in a quantum computer. Preliminary in- 
vestigations focused on the nature of the QWs themselves. For example, Kempe 
£D has shown that the hitting time of the discrete QW from one corner of an 
u( " A-bit hypercube to the opposite corner is polynomial in the number of steps, 

. n, whilst it is exponential in n the classical case. Subsequently Shenvi et al. jS] 

showed that a QW can perform the same tasks as Grover's search algorithm, 
and Childs et al. [S| introduced an algorithm for crossing a special graph expo- 
nentially faster that can be done with a classical RW. Kempe has recently 
reviewed the field. 

In the classical RW on the line, the "walker" (the particle or system per- 
forming the RW) randomly takes one step to the right or to the left depending 
on the result of tossing a coin. After n steps, the probability of finding the 
walker at a distance m from the origin is given by the binomial distribution, a 
Gaussian for large n with a standard deviation a — ^/n. In the QW, the role of 
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the coin is played by a qubit (as, e.g. a two-level atom or a spin-1/2 particle). 
As its classical counterpart, the quantum walker moves to the right or to the 
left depending on the internal state of the qubit. After each displacement, the 
state of the qubit is set to a superposition state by means of a suitable unitary 
transformation, typically a Hadamard, that plays the role of the toss of the coin 
in the RW. Yet the QW is not a random walk, as its time evolution is completely 
deterministic. The probability distribution in the quantum case is very different 
from the classical one: it resembles the Airy function (Fig.l) and has a standard 
deviation that is linear with n. This is the discrete time QW that should be 
distinguished from the continuous time QW [21E]> which we will not consider 
here. 

Possible implementations of the QW have been proposed by a number of 
authors 111! IT^] . Here we show that a classical implementation of the QW 
is possible, in analogy with other processes usually associated with quantum 
computing ^JEHEIE]- Indeed, we point out that a classical implementation 
very similar to the one we are proposing has actually been implemented by 
Bouwmeester et al. ! m the context of the optical Galton board, without the 
authors explicitly noting this. Other classical (interferometric) implementations 
of the QW have been proposed recently [HI E9, but in them the number of 
necessary optical elements grows quickly with the number of steps in the QW, 
something that does not occur in our scheme. Finally, by re-examining the 
difference equations for the walker we show that the nature of propagation is 
simpler than has been previously appreciated. 

In our classical approach the role of the walker is played by the frequency 
of a light field, and the role of the coin is played by its polarization state. The 
light field 

i 

E = E m cxp [i (u>o + mu>) t — ik m z] + c.c, (1) 

m=—l 

(loq is the carrier frequency and uj is the frequency difference between successive 
frequency components), can be represented by the abstract state 



IVO = 22 [Rm\m,x) + L m \m,y)}, (2) 

rn— — I 



where R m = x-u m and L m = y-u m (u m = E m / E m ) and J2 m =-i + \L m \ 

1; the "basis vectors" \m, c) label the frequency and polarization, with c = x, y; 
we associated x(y) with the coin head (tail). 

To implement the walk, we require a unitary operator that performs V 



' v 

m±l,yj. The operation V can be physically implemented, e.g., with an elec- 
trooptic modulator (EOM) to which a linearly time dependent voltage is applied 
in such a way that the x (y) polarization component of the field frequency com- 
ponent (uio + mid) will see its frequency increased (decreased) by an amount 
Q. 
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After each jump in the frequency of the field, a Hadamard transformation, 
H Tn, x ^j = H m ' x ) l m ' V)\ nas ^° ^ e implemented. This can be done opti- 
cally by means of a half-wave plate (HWP) with its fast axis forming an angle 
7r/8 with respect to the x axis ^3 El- Finally, the QW is implemented by 

the repeated action on the state of the operator HV, i.e. after n iterations 

r - ~ i 71 

— HV \ip (0)), that can be written 

+ N 

M n )) = [R-m,n\m,x) + L m . n \m,y)} (3) 

m=— JV 

Rm,n — {Rm-l,n-l + £m+l,n-l) , (4) 

Lm^n = — y= {Rm—l,n—l imtl,n- 1) 7 (5) 

where i? mj o = i m ,o = if m 7^ and R m .~i = i-m.-i = These are the 
standard QW equations. Finally, the intensity of each frequency component of 
the light field, which is the optical analog of the probability of finding the walker 
at position m at iteration (time) n, is given by P m .n = \Rm,n\ 2 + |^m,n| 2 j which 
is represented in Fig.l. 

In order to implement n steps the best option is to introduce the described 
elements in an optical cavity, Fig. 2. The cavity imposes a constraint that the 
optical frequencies must fit within its set of eigenfrequencies. Thus, the time 
dependent electric field applied to the EOM and the cavity length must be 
adjusted in such a way that the frequency shift Co = /lufsr with lofsr the 
cavity free spectral range and / an integer number. Consider, e.g., that a light 
pulse with a spectral width Alu is initially injected in the cavity. Then, in order 
to perform a step of the QW at each cavity roundtrip, the step size Co must be 
large enough to avoid significant overlap between the spectrum of the displaced 
pulses thus the frequency steps being well resolved. The total number of steps 
that can be made in this QW then depends on both cavity losses and EOM 
bandwidth. 

The experiment of Bouwmeester et al. |19j can be seen as a realization of 
the QW very similar to the one proposed here. These researchers proposed 
and studied, both theoretically and experimentally, an optical implementation 
of the Galton Board (the quincunx) . What they actually implement is a grid of 
Landau-Zener crossings through which a light beam propagates, and concentrate 
on the study of recurrences in the light spectrum. A simplified version of their 
experimental device is that represented in Fig. 2, but with the QWP replaced by 
a second EOM with its axis rotated ir/4 with respect to the first EOM, which 
introduces a dephasing between the two polarization components. Although this 
unitary operation does not correspond to a Hadamard transformation, it can 
be shown that it leads to an essentially identical QW (details to be reported 
elsewhere). The main difference with our proposal is that the frequency shift 
introduced by the EOM is smaller than lofsr and then each step in the QW 
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takes several cavity roundtrips. In Fig. 6 of their paper the QW is clearly 
seen. Bouwmeester et al. considered this case as a demonstration of the 
coherence quality of their system, and did not note its significance to QWs; 
their focus on the observation of recurrences in the spectrum led them to study 
other aspects of their system. 

Let us now re-examine the linear difference difference equations 111 !■")[) . They 
admit a formal solution that has been studied from a number of points of view, 
usually with a focus on identifying its asymptotic behavior for large n HO E] 
as it allows for the extraction of much information. Nevertheless the formal 
solutions presented to date do not rely explicitly on a crucial feature of Eqs. 
H4I5|) . which we now explain, that greatly simplifies a physical understanding of 
their solution. A little algebra reveals that the solutions R m ,n and L m>n of (|4I5|I 
also satisfy 

a m,n+l — a m,n-l + ~/= [ a >m— l,n — a m+l,n] j a = R,L. (6) 

This is a remarkable equation, since it demonstrates a dynamical independence 
of the evolution of the two coin states R and L. Thus there are two essentially 
independent walks, coupled only by the first step that links the a m ,i to the a m .o. 
After that the two walks can be studied independently of each other. 

The most naive continuous limit of © would involve a first derivative with 
respect to time and a first derivative with respect to space, and would suggest 
waves propagating only towards +00 for both R and L, in apparent violation 
of the symmetry of the problem. But this is too simplistic, given that for both 
R and L one can look for solutions of the form, 

a m ,n = An,ri + ( — 1) \n,n CO 

where satisfy 

^m,n+l — ^m,n-l = ^7/^ [^m-l,n — ^m+l,n] i ($) 

restoring the symmetry. Of course, there is not a unique specification of the A^ 
in terms of the fundamental a m , since 

O m ,l = Am^l — ^m,l- 

The specification of the a m ,o and a TOj i, which completely specifics the initial 
conditions required to solve JHJ), does not suffice to determine the initial condi- 
tions A^ o and A^ 1 required for the solution of (JHJ uniquely. Nonetheless, it is 
possible to rigorously develop the solution of the equations © in terms of the 
fields A^ n ; this we defer to a later publication. The point we wish to stress here 
is that in the limit of A^ n that are slowly varying in n and m we can introduce 
continuous functions A ± (x,t) and understand JSJ as the discretization of the 
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differential equation 

V ( m > d a± (r f y_ -r 1 y ( /Xx ) ° A ± (rt) (~\0) 

2^ (2k + l)\dt 2k + 1 [ ' j_T V2^ (2* + 1)1 9x 2fc +i 1 ' ; ' U j 

& — A.* — 

where At and Ax denote the temporal and spatial increments, respectively. 
Keeping only the first two terms and approximating the third derivative in time 
using the equation at the lowest order we obtain 



d 1 &' 



A ± (^r) (11) 



where r = t/ At and £ = x/Ax. In this slowly-varying approximation the con- 
ditions {HI can be approximated as 

A ± (m(Aa;),0) w A± fi = - (a m>0 ± o m ,i) , (12) 

where we have made use of Eq.0 and assumed that a mj i = 1 — A~ 1 w 
— A~ . Thus in this limit (|12l) provides the initial conditions for <|llfl 
and those equations can be solved by Fourier analysis. There are then two 
fields A^ (£, r) that can be associated with each side of the coin. This feature 
persists when the rigorous solution is constructed in this terminology, where 
there is a (temporal) "ferromagnetic" field A+ n and an "antiferromagnetic" 
field (—1)™ j4^ iB for each coin side. 

Returning to Ea. (|ll(l we make use of Eas. p2|l and take as initial conditions 

A ± (£, 0) = a , G (0)±a_ M G (-l)±a M G (1), with G (£„) = A^exp [- (£ - £ ) 2 / (2a) 2 

and AT a normalization factor; here we build in the fact that (jl 1|> is only correct 
for the long wavelength components by taking an initial condition that "smears 
out" the lower wavelength components. The solution is easily found analytically 
(23 E| for A^ (£, t), and we can write the final result (up to a normalization 
factor) for both R and L as a (£, r) = A+ (£, r) + (-1)" A~ (£, t) , with 

A ± (tr)=a , a Z(±^r) (13) 
± a-i tl Z (± (£ + So) - t) ± oi,i^ (± (£ - So) , r) , 

* = bits «P ( 3g2 J ^ j > ( 14 ) 

where A = t-r/y/2, B = t/(4V§), C = a 2 and A, (at) is the Airy function |2T]; 
the R and L solutions differ only through the different values of a m> o and a mj i 
appearing in 113( 1 . The appearance of Airy functions in the full solutions of (1 II oil 
|5] can thus be understood as associated with the form of equation Ijlll) , which 
workers in fiber optics will recognize as the classical equation for the propagation 
of light in a fiber with no group velocity dispersion but a third-order dispersion 
term. The linear dependence of the standard deviation on n arises, of course, 
simply because of this propagation. Solution lfH3|) is represented in Fig. 1(b) for 
a = 0.4, and the similarity with the QW in Fig. 1(a) is clearly apparent. 
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In conclusion: Wc have shown that the QW along a line can be simulated in a 
purely classical implementation, involving nothing more than wave interference 
of electromagnetic fields. And, indeed, it has in fact already been simulated in 
the laboratory in the work of Bouwmeester et al. ■ 

Further, this classical nature of the propagation is perhaps not surprising. 
After all, the standard QW is a generalization of the quantum mechanical prob- 
lem of a spinless particle with hopping amplitudes between sites, familiar from 
solid state physics if the time variable is continuous. That latter problem, which 
gives a simple Schrodinger equation in its continuum limit, is clearly classical in 
the nature of its propagation, as attested to by the appearance of the Schrodinger 
equation in classical beam propagation problems. The generalization involved in 
concocting the standard QW problem is the inclusion of a spin variable. What 
we have shown here is that this generalization does not affect the dynamics in 
an essential way. Except for an initial coupling in the first two time steps, the 
evolutions of the amplitudes associated with the two sides of the coin proceed 
independently. 
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Figure Captions 

Figure 1: (a) Probability distribution for n = 200 for both the classical 
(dashed) and quantum (continuous) random walks. The initial conditions cho- 
sen for calculating the QW were i?o,o = I/a/2 and Lo.o = i/V%, see Eas. (l4l5|l . 
Notice that the quantum P m is null for odd m at odd n. We have represented 
only nonzero values, (b) Continuous limit of the QW as given by Eq. ((T3|l for 
a = 0.4 and t = 200 with the same initial conditions as in (a). 

Figure 2: Scheme for the optical implementation of the QW in a Fabry- 
Perot cavity. The electrooptic modulator (EOM) shifts the field frequency up 
or down in u>/2 depending on its polarization and a quarter- wave plate (QWP) 
with its axis forming an angle tt/8 with respect to the a;-axis, performs the 
Hadamard transformation (notice that the light passes twice through each in- 
tracavity element every roundtrip). 
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